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We show that systems with negative specific heat can violate the zeroth law of thermodynamics. 
By both numerical simulations and by using exact expressions for free energy and microcanonical 
entropy it is shown that if two systems with the same intensive parameters but with negative 
specific heat are thermally coupled, they undergo a process in which the total entropy increases 
irreversibly. The final equilibrium is such that two phases appear, that is, the subsystems have 
different magnetizations and internal energies at temperatures which are equal in both systems, but 
that can be different from the initial temperature. 
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The zeroth law of thermodynamics is among the most 
fundamental assumptions concerning macroscopic sys- 
tems in equilibrium. In its textbook form, it concerns 
three systems A, B and C at equilibrium, that is, such 
that their macroscopic variables are time-independent: 
if A and B are in thermal equilibrium with each other, 
that is, if no heat flow arises when they are brought into 
thermal contact, and if further system B is similarly in 
thermal equilibrium with system C, then the zeroth law 
states that system A will also be in equilibrium with sys- 
tem C . By bringing two systems in thermal contact, one 
means allowing a weak coupling through which energy 
can be exchanged between the systems. From the above 
statement of the zeroth law follows the existence of an 
intensive variable, the temperature, which serves to pre- 
dict the behaviour of two initially isolated systems at 
equilibrium when they are brought into thermal contact. 
It states that the equality of the temperature parameter 
is a necessary and sufficient condition that no irreversible 
heat exchange will occur. 

It is of considerable interest to fully understand the 
possible limitations of such a basic law. In the following, 
we shall show that systems with negative specific heat 
can indeed violate the zeroth law. This is more than 
just an academic concern: systems with negative specific 
heat do exist in Nature. Although it is readily shown that 
systems well described by the canonical ensemble, such 
as those which are weakly coupled to an environment 
at fixed temperature, cannot have negative specific heat, 
others may. In particular, if a system is entirely isolated, 
that is, if it is in the microcanonical ensemble, the specific 
heat may well become negative To this end, how- 
ever, it is obviously necessary that the microcanonical 
ensemble lead to significantly different predictions from 
the canonical one. This cannot happen for large systems 
with short-range interactions. However, if either the in- 
teractions are long-range or, perhaps more realistically, if 
the system is small with respect to the interaction range, 



then it is not possible to refer to general statements con- 
cerning the equivalence of ensembles and negative specific 
heats can be, and indeed have been, observed. Such sys- 
tems include among others gravitational systems (long- 
range interactions) and atomic clusters (small systems) , 
see @, I, i, i, 1, U, B 0, M (U 0, H El . Analytically, 
mean-field models, involving interactions of all particles 
with all, are well-known to display negative specific heats 
0, ITEj in the microcanonical ensemble for parameters 
which, in the canonical ensemble, would correspond to 
the vicinity of a tricritical point [16| . 

Since systems with negative specific heat are thermo- 
dynamically unstable when they are thermally coupled 
to the surrounding medium [l7| , anomalous behaviour is 
surely to be expected when such systems interact. How- 
ever, it is not apparent that this will cause violations of 
the zeroth law. The reason is precisely that the zeroth 
law always involves coupling between systems. Heat ex- 
change is thus always allowed whenever the zeroth law is 
tested, and one might hence think that this always leads 
to a canonical-like case, in which the restriction on fixed 
energy is lifted and negative specific heats cannot occur. 
We shall show that this is not the case. 

In order to test the zeroth law for small systems in the 
microcanonical ensemble, one also faces yet another issue: 
for such systems it is often not straightforward how to 
define the temperature posited by the zeroth law. Under 
these circumstances, how, then, shall we know whether 
two systems really have the same temperature or not? 

A simple way around the last question is straightfor- 
ward: it suffices to take two identical systems, that is, 
two systems having the same Hamiltonian with the same 
energy, volume and number of particles, and then to cou- 
ple them weakly. Under these conditions, whatever the 
temperature is in each system, it certainly will be the 
same for both; yet irreversible changes may occur in the 
composite system even under such circumstances. 

To be specific, let us focus on the following Hamilto- 
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nian describing a system of N classical XY-rotors with 
phases defined by 
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where J > is the global ferromagnetic coupling and 
K is the nearest neighbors coupling; it can be negative 
or positive. The rotors are placed on a one-dimensional 
lattice with periodic boundary conditions. When K = 0, 
Eq. fl} is the classic Hamiltonian Mean Field model in- 
troduced in [l8j |. which has a second order phase tran- 
sition at the critical temperature T c = J/2 [l8| in the 
canonical ensemble (though the microcanonical and the 
canonical ensembles are equivalent when K = 0). This 
phase transition is characterized by the behavior of the 
magnetization order parameter m, 
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For K < it was shown [15[ that this system shows in- 
equivalence of ensembles for an interval of values of K. 
Specifically, there exists a region in which this system 
shows negative specific heat in the microcanonical en- 
semble. 

The canonical equilibrium properties are determined 
by the free energy per particle f(0) [l5( (J = 1), 
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where A(z, a) is the largest eigenvalue of the transfer ma- 
trix given by the following operator 
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In the microcanonical ensemble the thermodynamic pro- 
perties are given by the entropy per particle s M that is 
obtained from ([3]) by using the mean-field formalism in- 



troduced in Ref. [16] 
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where e is the energy per particle. The canonical entropy 
s c is found from l[5]) by a Legendre transformation with 




FIG. 1: Caloric curves (T vs e) in the canonical (black) and 
microcanonical (red) ensemble for the model JT) obtained from the 
exact expressions J3j and J5j, K = —0.178. 



respect to and has therefore the usual concavity prop- 
erty. The same does not hold, however, for s M (e) as is 
seen in figure [1] 

When we now couple both systems, the full system is 
still isolated and the total energy E = E\ + E2 + Ei nt is 
constant. As we consider weak couplings, the interaction 
energy can be neglected, so the total energy is the sum 
of the individual energies of each subsystem. Under this 
condition the total entropy per particle can be written as 



s(ei,e 2 ) = - Mei) + s M (e 2 )] 



(6) 



with s M given by ([5]) and the restriction e 2 = 2s— £1. Here 
e~ l is the energy density of subsystem 7 and s is the energy 
density of the full system, which is constant. The correct 
description of the full system is thus given by the mi- 
crocanonical entropy, although the energy conservation 
constraint was removed for each individual subsystems. 
By the second law of thermodynamics, the equilibrium 
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will be such that the total entropy 



will be maximum. So, the total entropy per particle is 
given by 
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From this optimization problem we obtain the following 
condition for the maximun entropy: T^(e\) — T 2 (2e — 
£*), where el is the energy that maximizes the total en- 
tropy and T^(e-y) is the temperature of the subsystem 
7 defined as the inverse derivative of the entropy with 
respect to the energy at constant N, V. 

By using ([5]) and |(7J) we can predict the equilibrium 
values of magnetization, energy density and temperature 
of each subsystem. In figure [2] we show s(e,ei) vs £1, 
microcanonical and canonical. For the microcanonical 
case there are two maxima. Because both subsystems 



FIG. 2: Total entropy per particle, s(e,ei), for the coupled sub- 
systems canonical (black) and microcanonical (red) , the dotted line 
indicates the value of energy density before coupling, for identical 
subsystems with negative specific heat. The values of the parame- 
ters are Ki = K 2 = -0.178, e = 0.55597. 



are identical, in principle which subsystem evolves to the 
higher energy is random. The total entropy before the 
coupling is now a minimum, that is, an unstable state. 
That this can occur was already noted in 0]. Therefore, 
the state of maximum entropy is now obtained when the 
subsystems have different values of magnetization and 
energy density, i.e. when two phases appear (a similar 
phenomenon was observed in [19(). The total entropy 
is therefore increased irreversibly if two subsystems with 
the same intensive parameters with negative specific heat 
are thermally coupled. Note that the final temperatures 
of both subsystems are identical, but in general differ 
from the initial temperature 20] . This difference is often 
too small to be detected numerically but its existence is 
rigorously established by the analytical results. 

To confirm our results we have performed numerical 
simulations of the complete system using a fourth-order 
symplectic algorithm with a time step 0.1 2l|. We run 
the simulations for a time interval r eq to let both sub- 
systes reach equilibrium without interaction (i.e. rj — 0). 
Once they have equilibrated, we increase the coupling lin- 
early during a time interval r a , after which the coupling 
is maintained constant, 77 > 0. 

We choose parameters in the region in which the un- 
coupled subsystems exhibit negative specific heat, and 
both subsystems have the same values of all their pa- 
rameters initially: energy per particle (si = £2), mag- 
netizations (mi = rri2) as well as particle numbers 
(JVi = A/2). From these equalities it follows that the 
derived variables, temperatures (T) and specific heats 
(c) are identical as well. The value of these ther- 
modynamic variables are defined through the micro- 
canonical entropy: T(e) — (ds^(e)/de)" 1 and c(e) = 



time 



— (d 2 s fl (e)/de 2 )~ 1 (ds fl (e)/de) 2 . The value of the magne- 
tization is such that the entropy for e fixed is maximum, 
therefore to = m(e) as well. Note further that for Hamil- 
tonians with a quadratic kinetic energy, as (|TJ) , we may 



FIG. 3: Temporal evolution of energy density (top) and magneti- 
zation (bottom) of subsystem 1 (black) and subsystem 2 (red), us- 
ing the coupling H? t . The curves were obtained by averaging over 
a sliding time window. The values of parameters are K\ = K 2 = 
-0.178, N 1= N 2 = 5000, -q = Q.l,N int = 10, e? = e§ 



sg = 0.55597. 



without appreciable error evaluate the temperature us- 
ing the expression for the kinetic temperature, which is 
given in this case as twice the mean kinetic energy per 
particle, since the finite N microcanonical corrections are 
negligible. This identification was verified by comparing 
with the momentum distribution, which we found to be 
Maxwellian. The measurement of the temperature in the 
simulation is therefore straightforward in this case. 

Here we report the effect of the following choice for the 
coupling Hamiltonian 



11 int ~ 
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where p] is the momentum of the i—th particle of subsys- 
tem 7. The advantage of this apparently peculiar choice 
is that this coupling does not exert a direct influence on 
the magnetizations of the two subsystems. 

In Figure [3] we show the behaviour of the internal en- 
ergy and the magnetization of both subsystems. Clearly, 
they do not remain at their initial equilibrium values; 
rather, due to finite size effects, the systems jump be- 
tween the two degenerate equilibrium states of the cou- 
pled system 22j: thus either subsystem one is in a mag- 
netized state and subsystem two in an to = state, or 
viceversa. The oscillation does indeed roughly take place 
between the two equilibrium values of the magnetization, 
which for the parameter values used here can be com- 
puted, along the lines described above, to be T* = T 2 = 
0.25, el = 0.56594, e* 2 = 0.546, mj = 0,m* 2 = 0.32. There 
is a problem, however, with the fact that neither subsys- 
tem ever comes reasonably close to being truly paramag- 
netic (to = 0)[13|. As it turns out, this is also a finite-size 
effect. As is readily seen, the mean-field Ising model at 
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criticality satisfies 
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which states that the typical magnetization fluctuations 
at, or in the vicinity of, a second order phase transition 
are of order iV -1 / 4 . A study of the TV-dependence of the 
magnetizations observed in our system is compatible with 
this scaling behavior [2(J. Thus, these large fluctuations 
due to the vicinity of a second order phase transition, 
which is known to exist at e c « 0.5633, are the reason 
why the zero magnetization phase is not clearly observed 
in our finite systems. 

The reason for the appearance of two phases lies, of 
course, in the thermodynamic instability of systems with 
negative specific heat, as discussed previously. Such sys- 
tems can only exist through the existence of constraints, 
such as the conservation of energy, which keeps them 
from relaxing to a more probable state. What we have 
therefore shown is simply that coupling two such sys- 
tems with each other provides enough freedom for an 
irreversible relaxation process to take place. It can also 
be shown that this phenomenon is not limited to 
the specific case studied here: any system with negative 
specific heat will show similar behaviour. The type of 
coupling used is also unimportant. 

Further, it is possible to use a third system (with posi- 
tive specific heat) as a thermometer to measure the tem- 
perature of the systems with negative specific heat. Thus, 
we can check that our systems are in a stable equilibrium 
with the thermometer at the same temperature before we 
put them in thermal contact. There exists, however, a 
condition over the third system: it needs to have a small 
heat capacity, otherwise it may allow large enough energy 
fluctuations to drive the system with negative specific 
heat out of the microcanonical ensemble (this condition 
on the thermometer was already discussed in 0, @, 0]; 
explicit simulations of this situation will be presented in 
[201 ] . where we use an ideal gas as the thermometer). 

Summarizing, we have shown how coupling two sys- 
tems with negative specific heat leads to an irreversible 
change in the intensive variables of the subsystems. In 
the simulations, the entire system does not settle down to 
a well- defined state. Rather, it displays slow oscillations 
between two degenerate states of thermodynamic equi- 
librium due to the finite size of the systems. It should 
further be noted that neither of these states corresponds 
to the canonical equilibrium state. 

Since both subsystems began at equilibrium with the 
same intensive parameters, the validity of the zeroth law 
implies that the coupling should not produce any no- 
ticeable effect. Since, as we have seen both through an 
analytical approach and by explicit numerical work, the 
system in fact relaxes to an inhomogeneous state in which 



the two weakly coupled subsystems have widely different 
values of the magnetization and of the internal energy, we 
have unambiguously observed a instance of a violation of 
the zeroth law of thermodynamics. 
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Negative specific heat may arise in systems that are of 
mean-field type as defined in fl(|, near the tricritical 
point. For this kind of systems we can write a Landau- 
like free energy density and we can approximate it by 
<p — a-zm 2 + aim 4 + a^m 6 . If the system is near the 
tricritical point (02 « 0, 04 ~ 0) then m ~ N~ a with 
a — 1/6, but if the system is near a critical point, 02 ~ 
and 04 is non zero, then a = 1/4. In either case, for very 
large N the normal a — 1/2 behaviour is recovered. 



